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Abstract. We discuss the 3-window modeling of cold, dense QCD matter equations of state at density 
relevant to neutron star properties. At low baryon density, ub < 2ns(n^: nuclear saturation density), 
we utilize purely hadronic equations of state that are constrained by empirical observations at density 
Ub ~ Us and neutron star radii. At high density, ub > 5ns, we use the percolated quark matter equations 
of state which must be very stiff to pass the two-solar mass constraints. The intermediate domain at 
2 < ub/tIs < 5 is described as neither purely hadronic nor percolated quark matter, and the equations of 
state are inferred by interpolating hadronic and percolated quark matter equations of state. Possible forms 
of the interpolation are severely restricted by the condition on the (square of) speed of sound, 0 < < 1. 

The characteristics of the 3-window equation of state are compared with those of conventional hybrid and 
self-bound quark matters. Using a schematic quark model for the percolated domain, it is argued that the 
two-solar mass constraint requires the model parameters to be as large as their vacuum values, indicating 
that the gluon dynamics remains strongly non-perturbative to ub ~ lOris. The hyperon puzzle is also 
briefly discussed in light of quark descriptions. 
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1 Introduction 

Neutron stars are the cosmic laboratories for the studies 
of quantum chromodynamics (QCD) at low temperature 
and high density [UEIE]. Solving Tolmann-Oppenheimer- 
Volkoff (TOV) equation together with the QCD equations 
of state, we can calculate the neutron star mass-radius 
(M-R) relation which can be confronted with the obser¬ 
vations. At given conditions such as rotation, magnetic 
field, temperature, each equation of state gives the unique 
M-R relation. This procedure is invertible; using Lind- 
blom’s algorithm |3], we can extract an equation of state 
uniquely from a given M-R curve [5]. 

For conventional neutron stars, the shapes of M-R 
curves are strongly correlated with pressures at several 
fiducial densities piT] . Pressure around ~ 1 —2ns (rb'- 
baryon densit 30 Ug ~ 0.16fm“^: nuclear saturation den¬ 
sity) determines the overall size of typical neutron stars, 
and beyond ~ 2ns the M-R curve moves toward the ver¬ 
tical direction without large variation of the radius. Its 
slope is correlated with pressure at n^ 2 — 4ns. The 
curve eventually ends at the maximal mass which is pri¬ 
marily determined by pressure at ns > 4ns. Thus the 
slope and maximum mass carry most direct information 


^ Notations: [ub, fj,B,-'') with subscript B are used for 
baryonic quantities, and (n, /r, ■ • •) for quarks. 


for dense matter beyond the conventional nuclear matter 
domain. 

In this respect, the discoveries of two-solar mass {2Mq) 
neutron stars [Slin] have significant impacts. The existence 
of such stars mean that the QCD equations of state at 
Ub > ‘2ns are very stiff, i.e., the pressure (P) at given 
energy density (e) is very large to resist the gravitational 
collapse. In fact a variety of soft high density equations of 
state have been ruled out by this mass constraint alone. 
The neutron star radii should put further constraint nni 
HUliailslIIllllSlIle], although the determination contains 
more uncertainties than the mass measurement. 

There are several purely nucleonic equations of state 
that can pass the 2Mq constraint. But typical hadronic 
models, in their typical versions, indicate that hyperons 
are likely to appear at ns > 2 — 3ns, and they greatly 
soften the equations of state; in general the massive de¬ 
grees of freedom, near its mass threshold, adds mass times 
its number density to the energy density, but adds only 
little amount of pressure. The conflict between the 2Mq 
constraint and hyperon softening is called ’’hyperon puz¬ 
zle” [T7iiT^rm[^iTr| . 

Similar softening problems occurs for typical hybrid 
equations of state, in which hadronic and weakly inter¬ 
acting quark matter are treated as distinct and the full 
equation of state is made through the Maxwell construc¬ 
tion for the first order phase transition. Typical models 
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■ Few quark exchange ■ many-quark exchange ■ Baryon overlap 

■ structural change of hadrons ' Quark Fermi sea 



Fig. 1: The 3-window description for QCD matter, (i) At ub < 2ns, nucleons are well separated and exchange only few quarks 
(or mesons), (ii) At ns > 2ns, baryons other than nucleons start to appear and baryons exchange many quarks, developing 
many-body forces. Accordingly there are structural changes in hadronic wavefunctions. (iii) At ns > 5ns the baryons begin to 
percolate and the quark Fermi sea develops. Until ns reaches ~ lOOus, the matter remains in the strongly correlated regime. 


predict the phase transition occur at ns ~ 2 — lOug. Since 
the first order transition has a jump in energy density at 
fixed pressure, it always softens the equations of state. 
This calls debates on the existence of quark matter at the 
neutron star core [SUES]. 

We shall take up these problems as great opportuni¬ 
ties to study a matter which is neither purely hadronic 
nor weakly interacting quark matter. To be clear, we em¬ 
phasize that the purpose of this paper is not to predict 
the equations of state. Rather, our primary purpose is to 
delineate the properties of strongly correlated QCD mat¬ 
ter and develop a plausible picture for non-perturbative 
QCD at finite density. This will be done by examining 
supposed equations of state that are being more and more 
constrained by the neutron star observations and the QCD 
predictions from the chiral effective theory (ChEFT) and 
perturbative QCD (pQCD). Then we will argue which as¬ 
pects of the QCD matter are relevant to resolve the afore¬ 
mentioned softening problems. 

To infer the properties of such matter, we first note 
that the softening problems occur somewhat outside of the 
domain of applicability of the conventional descriptions. 
In this respect, we should refer to the chiral effective the¬ 
ory [^E5lE5ll27[l28lE^ and pQCD [5011^15^1^ . These 
methods can provide not only quantitative predictions but 
also predict the domain of applicability; in particular these 
methods can address when the underlying pictures behind 
the calculations require modifications. The hyperon prob¬ 
lem occurs around ~ 2 — within purely hadronic 
descriptions. But such descriptions, at ns > 2ns, may 
cause problems, or at least raise questions, on the con¬ 
vergence of many-body forces, appearance of new bary- 
onic degrees of freedom, and possible structural changes 
in the QCD vacuum and hadronic wavefunctions. To un¬ 
derstand all these issues, it is important to go back to more 
microscopic description from which effective theories are 
derived. On the other hand, the quark matter problem oc¬ 
curs at Ub ^ 2 — lOris within the description of weakly 
interacting quark matter. But at such density the weakly 
interacting picture for quark matter is not quite valid as 
can be seen from pQCD equations of state; the current 
state-of-art calculations to 0{al) clarify that the inter¬ 


actions are important already at relatively large quark 
chemical potential p, ~ 1 GeV or ns ~ 10^77,^. Moreover 
the pQCD at ns < IQ^Ug results start to show fairly 
large renormalization scale dependency indicating that 
in quark descriptions the matter at ns — 2 — lOno should 
be regarded as strongly correlated matter. 

These observations bring us to the 3-window descrip¬ 
tion for the non-perturbative domain of QCD matteiQ 
which was introduced by Masuda, Hatsuda, and Takat- 
suka [3S], and was also discussed in Ref. in 

slightly different versions. Below we will follow the descrip¬ 
tions of Ref. |38j , whose schematic picture is summarized 
in Figig We use purely hadronic equations of state at 
low density where baryons exchange only few mesons (or 
quarks). To be specific, at ub < 2ns, we use the Akmal- 
Pandharipande-Ravenhall (APR) equation of state, which 
is purely nucleonic [39]. At densiW high enough for the 
overlap of baryons or percolatiorjj [ID] , we construct the 
equations of state using effective quark models at the QCD 
scale Aqcd — 0.2 GeV. To be specific, at ub > 5ns, we use 
the 3-flavor Nambu-Jona-Lasinio (NJL) model plus vector 
and diquark interactions. Consequently we are left with 
the domain of 2ns < tt-b < 5ns for which appropriate pic¬ 
tures are not known even at qualitative level. Such matter 
is supposedly described as hadronic matter about to per¬ 
colate, or quark matter about to be confined. We infer the 
equations of state by interpolating the APR and NJL pres¬ 
sures as functions of chemical potential. In this work we 
will use smooth interpolation assuming the hadron-quark 

^ Or the convergence problems at p < 1 GeV can be clearly 
seen from the studies of Ref. [30] where the authors compared 
the 0 ( 03 ) and 0{a^) corrections. 

^ More modest interpolation approach is to connect hadronic 
equations of state to pQCD as worked out in Refs. [D^I33II34| . 
While these studies have more predictive characters, our pri¬ 
mary concern is how to interpret such pressure curves. These 
two approaches are complementary. 

If we take the volume of a baryon to be V = A-kTs/S 
with the radius r^, the baryons overlap at ns — 3 — 12n3 for 
Ts ~ 0.5 — 0.8 fm, assuming them to be static objects. The per¬ 
colation would occur at lower density because their momenta 
are not quite small. 
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crossover although the 3-window construc¬ 

tion itself does not reject the possibility of the first order 
phase transition as far as the jump of the energy density 
is not too large. For simplicity, we use a polynomial form 
for the interpolated pressure and require it to smoothly 
match with the APR and NJL pressures up to the second 
order of derivatives. 


In addition to the boundary conditions, the interpo¬ 
lated equations of state must obey several physical condi¬ 
tions. The interpolated pressure must (i) be stiff enough 
to allow the existence of 2Mq neutron stars; (ii) yield the 
square of speed of sound, = dP/de, which is positive 
and smaller than the square of speed of light, 0 < < 1. 

The reasonable interpolation satisfying these conditions is 
possible only for a restricted class of quark matter equa¬ 
tions of state at high density. We will study what kind 
of interactions are relevant for such quark matter equa¬ 
tions of state, and will show that such interactions follow 
from the knowledge of hadron spectroscopy and nuclear 
physics, provided that the gluon dynamics at finite quark 
density is not significantly different from that in the QCD 
vacuum 


This paper is structuralized as follows. In Secwe be¬ 
gin with preliminary remarks. We present graphical anal¬ 
yses of P(/i) curves which allow us to get rough but quick 
insights on the stiffness and speed of sound for given equa¬ 
tions of state. This method appears to be very useful to 
classify several types of equations of state. In Sec|^ we 
examine the characteristics of the 3-window equations of 
state and argue the quark matter softening problem. We 
discuss why the 3-window construction can give a stiff 
high density equation of state, by allowing quark matter 
equations of state stiffer than those typically appearing in 
the conventional hybrid construction. The observation is 
especially important when we consider soft hadronic equa¬ 
tions of state at low density. The self-bound quark matter 
equations of state are also discussed for further compar¬ 
isons. In Sec|^ we briefly explain a model of nuclear and 
crust equations of state. In Secj^we introduce a schematic 
quark model for percolated quark matter. The model con¬ 
sists of the standard NJL model plus additional vector and 
diquark interactions for baryon and nuclear phenomenol¬ 
ogy. We examine the impacts of the additional interac¬ 
tions and show that these contributions are important to 
express the sensible equations of state. In Secj^ we con¬ 
struct 3-window equations of state and compute the M- 
R relations by solving TOV-equation. We will show how 
severely the conditions of speed of sound constrain inter¬ 
polated equations of state and thereby the range of model 
parameters. Secis devoted to conjectures and interpre¬ 
tations of our equations of state. We discuss the hyperon 
problem from the viewpoint of quark descriptions. Physics 
behind our stiff quark matter equations of state is also dis¬ 
cussed in the context of gluon dynamics at finite density. 
The Secj^is devoted to summary and outlook. 




Fig. 3: Stiff (Pi) v.s. soft (Pq) equations of state. At given 
pressure P*, the curve Pi gives less energy density than 
the Pq curve. 

2 Preliminary remarks 

Following Ref.[3H], we consider equations of state in the 
grand canonical ensemble and work with P(/i). In this 
choice of thermodynamic variables, the problems in quark 
matter softening will become most manifest. The thermo¬ 
dynamic ground state has the maximal pressure at given 
/r. If there are several candidates of phases, we can just 
plot them in the plane of P v.s. ^ and can easily iden¬ 
tify the phase transition point. On the other hand, the 
discussions of the stiffness (P v.s. e) require a little more 
efforts; we need to calculate n{^) = dP/d^ and then use 
the relation £(/r) = /in — P to extract e. 

In this section we discuss a graphical method which 
can be used to get quick estimate of stiffness from a P(/i) 
curve. The method is also useful to get rough insights for 
the conditions of the speed of sound. Thus using a P(/i) 
curve alone, we can see the relation among phase transi¬ 
tions, stiffness, speed of sound, and then can understand 
how these conditions put constraints one another. 


2.1 Stiffness from P(/i) curve 

Suppose some P(/i) curve. Our task is to estimate e at 
some value of chemical potential, /t*. To find out e* = 
£(/i*), we notice that the tangential line of the P(/t) curve 
at /i = /t*, with the slope n* = 9P/i9/i|^=^,, intercepts 
with the P-axis at e*. This can be easily seen from FigH 
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Fig. 4: Stiff (P 2 ) v.s. soft (Pq) equations of state. P 2 is a 
clockwise rotated curve of Pq- 

To see the utility, let us apply the method to several 
examples. The first example is the bag model with the bag 
constant B{> 0). This constant must be added to (sub¬ 
tracted from) energy density (pressure) when the ground 
state becomes perturbative, because the perturbative vac¬ 
uum has larger energy than the non-petrubative vacuum. 
The pressure in the perturbative domain is given by P = 
co/r^ — B, and from which n = 4co/i^ and £ = Sco/i"*^ -I- B 
follow. Since P = {e — 4P)/3, larger B leads to softer 
equations of state. In the graphical method, this conclu¬ 
sion follows immediately as seen from Fig|^; at fixed value 
of P* = P(^*), a curve Pq, which has a larger bag constant 
B than Pi, leads to the larger length for the intercept for 
the tangential line. This means that Pi is stiffer than Pq. 

For another example, suppose a curve P 2 which is ob¬ 
tained by rotating Pq clockwise around the point where 
Pq = 0 (Fig 4). Clearly, the length of the intercept of 
the tangentialline gets smaller after the rotation, mean¬ 
ing that P 2 is stiffer than Pq. This sort of behaviors oc¬ 
curs when we include repulsive density-density interac¬ 
tions; more chemical potential (or energy density) is re¬ 
quired to increase the number density and pressure. The 
resulting P(/j.) curve grows slowly as /r increases. 

Later we will examine the impacts of various inter¬ 
actions on P(/i) curves. The results may be understood 
using the above two examples. 


2.2 Physical speed of sound from P(/r) 

Equations of state are under the constraints from the gen¬ 
eral relativity and thermodynamics; the square of speed 
of sound = dP/de must be positive and less than the 
square of the speed of light, 0 < < 1. We examine how 

to see these constraints in the P(/r) curves. 

To do so, first we notice the thermodynamic relations, 
dP = nd^j, and d£ = /rdn, from which one gets 

2 dP 91n/r 
® de 9 In n 


Fig. 5: The schematic plot for P(/r) curves with fixed 
sound velocities. At some chemical potential (indicated 
by the dot), we switch from some curve to fixed sound 
velocity curves. For equations of state to be causal, the 
P(/i) curves must grow faster than quadratic powers. 



Fig. 6: A unphysical P(/r) with an inflection point. The 
region with d'^P/dy? < 0 has the unphysical speed of 
sound, < 0. 


differential equation to get 


n = 



( 2 ) 


For instance, n oc. p? for an ideal gas with = 1/3, or 
n oc /r for the causal limit = 1, or n = const, for infinite 
speed of sound, = 00 . Note also that if n is a decreasing 
function of /r, the square of the speed of sound is negative 
and unphysical, < 0. 

In this example, the pressure at /i ~ p,* is (P* = P(/x*)) 
P ~ P, -I- const. X . (3) 

For an equation of state to be causal (c^ < 1), the pressure 
must grow faster than at any points on the P(/i) curve 
(Figjs]). 

In Ref. [38] . the authors imposed the positive suscepti¬ 
bility constraint, d'^P/dy? = dn/dfj, > 0. This is equiva¬ 
lent to requiring the speed of sound to be positive, because 


Suppose a small region around (n*,/j,*) where Cs is ap¬ 
proximately constant. In such a region, we can solve the 


d^P _ dn 
5^2 i9/r 


const. , i/cj-i 


( 4 ) 
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Fig. 7: A hadronic pressure Ph and quark equations of 
state, Pqo, Pqi, and Pq^- In the conventional hybrid con¬ 
struction, only Pqo can be used; the other curves are stiff 
(see Figsj^and]^, but cannot be combined with Ph- 

becomes negative only if < 0. The positive susceptibil¬ 
ity condition can be easily examined by checking whether 
P(p,) curves contain inflection points or nolQ (Fig|^. 

As we will see, the interpolation in the 3-window con¬ 
struction is tightly constrained by these conditions. In 
turn, the quark matter equations of state, which supply 
the boundary condition for the interpolated curves, are 
also severely restricted. 


3 Classification of quark equations of state 

In this section we classify the equations of state that con¬ 
tain quark matter at the neutron star core, and review how 
they pass the 2Mq constraint in different ways. The con¬ 
ventional hybrid, 3-window, and self-bound quark matter 
equations of state are discussed. 

3.1 The conventional hybrid equations of state 

In the conventional construction of hybrid equations of 
state, one uses hadronic and quark matter equations of 
state as distinct. At some density the matter undergoes 
the first order phase transition from hadronic to quark 
matter phase. To describe such situation, we must use a 
quark matter equation of state such that its pressure P(/r) 
approaches the hadronic one from below; quark pressure 
must grow sufficiently fast, otherwise we would not get 
the intersection point (Figj^. As we saw in the last sec¬ 
tion, such quark equations of state are typically soft, hav¬ 
ing troubles with the 2Mq constraint. This tendency has 
called debates on whether quark matter can exist at the 
core of neutron stars mm- 

There are several ways to bypass this stiffness prob¬ 
lem. For instance, we can allow stiff quark equations of 
state if hadronic equations of state extrapolated to high 
density are stiff [47ll48ll4^ . In that case the hadronic pres¬ 
sure grows so slowly that even a stiff quark pressure curve 

® To see this, in most plots we will not divide P(/i) by or 
the ideal gas pressure. 


Fig. 8: An example of stiff hybrid equations of state. As¬ 
suming the hadronic equation of state is stiffened at high 
density, stiff quark equation of state can find the intersec¬ 
tion point at low density. 



Fig. 9: An example of stiff hybrid equations of state. As¬ 
suming the intersection of Ph and Pq occur at very high 
density, hadronic equations of state, which are assumed to 
be stiff at very high density (like purely nucleonic equa¬ 
tions of state), build up most of the neutron star mass 
before the quark matter softening occurs. 

(which does not grow very fast) can have the intersection 
with the hadronic one (Figj^. In this stiff-stiff combina¬ 
tion, it is possible to pass the 2Mq constraint with large 
amount of quark core (although some cares are necessary 
for the causality condition). Another way to pass the 2Mq 
constraint is to use stiff high density hadronic equations of 
state and to locate the first order phase transition point at 
very high density [39]. In this case, the stiff hadronic mat¬ 
ter builds most of the neutron star mass before the quark 
matter softening occurs (Figj^. In this case the amount 
of quark matter core is tiny. 

While the above treatments can satisfy the 2Mq con¬ 
straint, there is fundamental weakness in the conventional 
hybrid construction. The problem is that we must com¬ 
pare two equations of state even though their domains of 
applicability might not overlap; purely hadronic picture is 
valid at low density while at high density its applicability 
is questionable; quark matter picture is valid at high den¬ 
sity but at low density it would be affected by large non- 
perturbative corrections such as confining effects. It is po¬ 
tentially dangerous to use high density hadronic equations 
of state to select out acceptable quark matter equations 
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Fig. 10: An example of 3-window equations of state. The 
boundaries are given at ~ 2ns and ns ~ for 
hadronic and percolated quark pressure, respectively. The 
3-window construction accepts pressure curves forbidden 
in the conventional construction, and allows us to explore 
a new class of (stiff) quark matter equations of state. 

of state. Also one cannot rule out quark matter equations 
of state even if their low density behaviors are at odd 
with nuclear phenomenology; the origin of problem might 
be just the wrong extrapolation of correct high density 
quark equations of state. 

Finally, let us note that the conventional construc¬ 
tion implicitly omits an equation of state which is nei¬ 
ther purely hadronic nor quark matter. The 3-window con¬ 
struction, which will be discussed next, will deal with this 
omitted possibility. 


3.2 The 3-window equations of state 

The 3-window construction is based on the picture that (at 
ns < lO^ns) QCD matter may consist of three domains; 
purely hadronic matter at low density, percolated quark 
matter at high density, and matter intermediate between 
these two matters. 

In practice, at low density we use equations of state 
purely based on the hadronic degrees of freedom. Hadrons 
are proper effective degrees of freedom to describe the sys¬ 
tem as far as the system is dilute enough. Baryons are 
well-separated and exchange only few mesons (quarks), 
so that the interactions are tractable and modifications 
of hadronic wavefunctions are negligible. Also there are 
only few excited baryonic states appearing in the virtual 
processes, and they can be manifestly taken into account 
or can be integrated out into a few numbers of effective 
interactions. These treatments, however, are difficult to 
justify at high density. In fact there are no established 
hadronic equations of state beyond ~ 2ns. Thus we will 
stop using the purely hadronic descriptions for density 
above ub — ‘2.ns. 

At high density, we use equations of state based on 
the percolated quark matter picture. Baryons overlap and 
quarks no longer belong to specific baryons. Then quarks 
develop the Fermi sea. The percolated matter picture, 
however, cannot be justified at density smaller than ns ~ 



Fig. II: The schematic plot for the equation of state in 
finite temperature QCD. By interpolating the hadron res¬ 
onance gas (HRG) and (resummed) pQCD pressure, the 
3-window pressure gives good qualitative description for 
the lattice data. 


5ns; three quarks gather to form a baryon. The confining 
effects become crucial so we will stop using the percolated 
quark matter descriptions for density below ~ 5ns. 

There is an intermediate region for 2ns ^ ns ^ 5ns 
where neither of the above pictures are applicable. The 
matter may be described as strongly interacting hadronic 
matter, or hadronic matter about to percolate, or quark 
matter about to be confined. To infer the corresponding 
equations of state, we simply interpolate the equations of 
state of purely hadronic and percolated quark matters. 

In the 3-window construction, we have more chances 
to get stiff quark matter equations of state than in the 
conventional hybrid one (Fig |10[ ). For instance, the extrap¬ 
olated hadronic equations of state at high density will not 
be used to rule out some of stiff quark matter equations 
of state. Also we will not judge the validity of a quark 
equation of state by its low density behavior because it 
should acquire huge corrections at low density. Any quark 
equations of state are acceptable as far as their behav¬ 
ior at high density is appropriate and provide the proper 
boundary conditions that allow the interpolated equations 
of state to be physical. These observations enable us to 
study a new class of quark equations of state which have 
not been fully explored, and the resulting quark equations 
of state can be stiff. 

In order to make the 3-window arguments more con¬ 
vincing, here we quote the finite temperature QCD equa¬ 
tion of state to use some analogy (Figj^ See also Refs. m 
[5T]j. At low temperature, the hadron resonance gas (HRG) 
model, which is a gas of non-interacting hadrons, gives 
very good description for the pressure of the lattice QCD 
data |521l531l541l55j . Around the (pseudo-)critical tempera¬ 
ture Tc for chiral (or deconfinement) transition, however, 
the HRG pressure starts to overshoot the lattice data. This 
is because around Tc thermally excited hadrons overlap 
and the interactions are not negligible; the inclusion of the 
interactions tempers the overshooting behavior. On the 
other hand, at high temperature the (resummed)pQCD 
theories give reasonable descriptions for the lattice data 
[M1I57] . Below 2 — 3Tc, however, the pQCD pressure tends 
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Hadronic 

(free) 


many-body forces 



Fig. 12: Our expectation for the equation of state in finite 
density QCD. By interpolating the hadronic and perco¬ 
lated quark matter pressures, we expect that the 3-window 
pressure provides reasonably good description. 


to overshoot the lattice datsj^ This is not surprising, be¬ 
cause the pQCD calculations do not include the confining 
effects which trap quarks into hadrons. The inclusion of 
the confining effects such as the Polyakov loop can tame 
the problem by suppressing artificial excess of the pQCD 
pressure [55] . 

Now notice that if we trusted the extrapolated HRG 
or pQCD pressure as real, we would construct the hy¬ 
brid equations of state in which we have pQCD gas at 
low temperature while HRC at high temperature! This 
senseless result tells us that one must be very careful 
when comparing the extrapolated equations of state. In 
contrast, if we stop using the HRC and pQCD gas pic¬ 
tures beyond the domains of their applicability, the do¬ 
main Tc < T < 2 — 3Tc is left over. We can infer the 
pressure in such domain by interpolating the HRC and 
pQCD gas pressure, and then get the reasonable result; 
the 3-window construction for finite T QCD equations of 
state works well. The domain Tc < T < 2 — 3Tc may be 
described as either strongly interacting HRC, or ”semi- 
QCP” [53], in which the quark-gluon plasma is subject 
to large non-perturbative effects. The finite density ver¬ 
sion of the semi-QCP is what we try to describe in the 
interpolated equations of state (Fig 12). 


3.3 The self-bound quark matter equations of state 

The thermodynamic ground state of a self-bound matter 
has nonzero number density 0 at threshold 

chemical potential /ith such that P(fith) = 0. The self¬ 
bound quark matter hypothesis states that the absolute 
ground state of QCD matter is (percolated) strange quark 
matter [60]. If so, at low chemical potential the strange 

® The overshooting behavior is clearly seen in pure Yang- 
Mills for T < 2 — ~ 0.5 — 0.8 GeV, while with quarks 

this tendency depends on the choice of the renormalization 
scale; the hard thermal loop results for LO, NLO pressure, and 
NNLO result with relatively small show the overshooting 
behaviors, but NNLO with large as gives the undershooting 
behaviors |57|. But the latter has large deviation from NLO 
pressure, and its validity is a bit uncertain. 


quark matter pressure is larger than the nucleonic one (see 
Pqi in FigQ and nuclear matter is just a meta-stable 
state. The decay from nuclear to strange quark matter 
is suppressed because of vastly different forms of these 
matters; wavefunctions do not quite overlap and it is hard 
to tunnel from one to another macroscopic state. 

Strange quark matter equations of state can be made 
very stiff, if /uth and n = dP /are small. In fact there 
are several examples of this sort that can pass the 2Mq 
constraint. 

However, there are important qualifications, if one tries 
to use strange quark matter to explain the 2Mq neutron 
stars. As we saw in the last section, to get stiff equations 
of state we wish to make pressure slowly growing function 
of /r. This condition requires n = to be small. 

However, if n is too small or quark matter is too dilute, the 
confining effects are not negligible. The confining effects 
turn dilute strange quark matter into hypernuclear matter 
by packing three quarks into a nucleon or a hyperon. 

To be more concrete, we consider the bag model, 

NcNi 4 

where we ignore the strange quark mass and as correc¬ 
tions. For this simplest version of the bag model, the star 
mass is [^ 


Afmax = 2.033 X 


/56MeV/fm3 

I B 


1/2 

Mq . 


( 6 ) 


So if i? > 56MeV • fm“^, the self-bound quark matter can 
pass the 2Mq constraint. In this case p,th = 274 MeV. But 
a question arises for the corresponding baryon density. It 
is 


nB(Aith) = 0.273fm-3 ~ I.7n, , (7) 

where the system is presumably too dilute to be free from 
the confining effects, with the fact that at this density 
purely nucleonic calculations do not find the serious prob¬ 
lems yet. There are also analyses with aj. corrections to 
the bag model, but the equation of state compatible with 
the 2Mq constraint has the similar onset density |5|. 

For these reasons we are inclined to think of self-bound 
quark matter as a candidate of small, light stars. In this 
case, we can take the onset baryon density at /ith to be 
high enough to apply the percolation picture. 


4 A model of hadronic matter: APR 

For the 3-window construction it is necessary to choose a 
hadronic equation of state. For very dilute region at ns < 
0.5ns, we use the Skyrme-Lyon (SLy) equation of state for 
the crust m- The crust part shares only small fraction of 
the neutron star mass, so the details are not essential for 
our discussions. On the other hand, at ng > 2ns there are 
several issues on the validity of existing hadronic equations 
of state, so we will not try to directly describe such region. 
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Then the remaining task is to choose one of equations of 
state for O.Sng < ub < 2ns, from the list available in the 
literature. 

For 0.5ns < ns < 2ns, the purely nucleonic equations 
of state are viable candidates. There are several nucleonic 
equations of state in the literature. They are constructed 
to reproduce the experimental data for symmetric nuclear 
matter at n^ ~ ns, but they are extrapolated to ns ^ 2ns 
in different ways. Also, the differences can be already seen 
in the extrapolation to asymmetric nuclear matter, and 
they affect the neutron star radii by ~ 1 — 3 km [6]. 

In this work we use the Avl8+(5t!+UIX* version of 
APR equation of state [39], that is based on variational 
many-body calculations for the non-relativistic potential 
models. The effective 2- and 3-body potentials are fixed 
by two nucleon scattering data below pion production 
threshold, spectra of light nuclei, and properties of sym¬ 
metric nuclear matter at saturation density. (Nowadays 
more QCD-based derivation of such potentials is avail¬ 
able in a systematic power counting in ChEFT . 

Then those potentials are combined with the advanced 
many-body calculations; for recent progress, see Refs. dzi 

I1H1I2SI-) 

The APR equation of state contains the first order 
phase transition for any proton fraction. The phase tran¬ 
sition is interpreted as the pion condensation through the 
two-body correlation functions. For the /3-equilibrium mat¬ 
ter, the transition occurs at ub — 1.4ns. (For symmet¬ 
ric nuclear matter, the transition density is higher, ub — 
2ns.) Although the strength of the phase transition is not 
too signihcant, we will see that this small amount of soft¬ 
ening already introduces difficulties on the construction of 
stiff equations of state at high density. 

Typically, the non-relativistic potential models provide 
softer equations of state at ub = 1 — 2ns than mean field 
type calculations. The mean held type hadronic equations 
of state are rather stiff at low density and can be com¬ 
bined with stiff quark equations of state at high density 
(see the discussions around Figj^. In fact there are sev¬ 
eral hybrid equations of state of this type which pass the 
2Mq constraint [4711481149] . On the other hand, with the 
non-relativistic model pressure at low density, the con¬ 
struction of stiff equations of state becomes more difficult 
problem [62] , especially when we use the conventional hy¬ 
brid construction. 

Several nucleonic equations of state, which are con¬ 
strained by the properties of symmetric matter at satura¬ 
tion density, can show considerable differences when they 
are extrapolated to asymmetric matter region. It is estab¬ 
lished that the neutron star radii are strongly correlated 
with the pressure at ns = 1 —2ns which is sensitive to the 
symmetry energy and its density dependence [SJUS]. For 
mean field models the radii are typically i? = 13 — 16 km, 
while in the potential model i? = 10 — 13 km. The recent 
estimates of the neutron star radii from the observational 
data seem to favor rather small radii of 9 — 13 km [lOl 
[IIl[IIl[I3[Il[ig with an exception suggesting R > 14 km 
m, although all these analyses contain some model as¬ 
sumptions. If these small radius estimates remain true. 


then we have to construct equations of state which are 
soft at Ub = 1 — 2ns and becomes stiff at n^ > 2ns. The 
soft equations of state at ns = 1 — 3ns also seems to be 
favored by the heavy ion collision data by Danielewicz et 

al. [55] . 

While the nucleonic equations of state are safely appli¬ 
cable only up to Ub ~ 2ns, the behaviors at low density 
already put severe constraints on high density equations 
of state through the boundary conditions, thermodynamic 
and causality constraints. Using the APR at low density, 
we need to combine the soft low density hadronic equa¬ 
tions of state with stiff dense equations of state. Such com¬ 
binations make the 2Mq constraints even severer. In the 
following sections we try to construct such equations of 
state. 


5 Modeling percolated quark matter 

In this section we will introduce a schematic model for 
quark matter descriptions, and examine the impact of var¬ 
ious effective interactions on the equations of state. 

5.1 Quark model 

At Ub < lO^ris, pQCD is not safely applicable and we 
must use some effective model for non-perturbative QCD. 
Our guideline toward proper effective models is hadron 
spectroscopy and nuclear physics. We consider the follow¬ 
ing schematic HamiltoniarQ 

n - Hnjl + + Hbb + • • •. (8) 

Let us clarify physics which we wish to express. 

"HNiL is the standard NJL model which is responsible 
for the chiral symmetry breaking/restoration and struc¬ 
tural change of the Dirac sea and quark bases. It is given 
by 

33njl = g(i7o7 ■ d + m- fijo)q 

G ^ 

-+ ( 917579 )^] + 8K{det{qRqL + h.c.) , 

(9) 

where m = diag.(TO„, m^, rus) is the current quark mass 
matrix, and tq = (2/A^f)^/^ and Ti(i = 1, • • • ,8) are the 
Gell-Mann matrices for flavors. The last term is Kobayashi- 
Maskawa-’tHooft interaction for the 174(1) anomaly [55] . 

^conf ^ expresses the confining effects that trap 3-quarks 
into a baryon. This term will be essential if we try to de¬ 
scribe hadronic matter microscopically. But right now we 
are not aware of the proper description in the practical 

^ The Hamiltonian in this paper is a bit simpler than that of 
Ref. [3H] in which the anomaly coupling between chiral and di¬ 
quark condensates is also included. The effect of such coupling 
can be mostly absorbed by changing the value of {Gv, H), so 
we will stick to the simpler version of the model. 
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models. Thus we will restrict the use of the model to the 
percolated domain, and then ignore assuming its 

contribution to be small. 

TLoge is a model for the one-gluon exchange contri¬ 
bution which is relatively short-range. We focus on the 
color-magnetic component: 

'^0GE = -H X! {q^C\-I^TA\A'q) 

A,A'- 2 ,b ,7 

+ i^TAXA'C'f') {q^CTA\A'q) ] , (10) 

where are the Gell-Mann matrices for colors. Here we 
kept only the attractive channel: color-flavor-spin anti¬ 
symmetric part. In hadron spectroscopy, it is responsible 
for the N — A splitting in the context of the constituent 
quark models. This interaction is called diquark-diquark 
interactions in researches of the color-superconductivity. 

"Hbb is inspired from baryon-baryon interactions. This 
is responsible for the quark version of the Walecka model 
[55] type descriptions. Since we have already included the 
cr-meson part in the NJL model, we add only the ui and (j) 
meson-like contributions or hard core type repulsive con¬ 
tributions 1551: 


S-\k) = 


(16) 


With these dehnitions, the diquark condensates (di, ^ 2 , da) 
correspond to (ds, sm, ud) quark pairings, respectively. 

The thermodynamic potential may be computed from 
the mean field particle propagators in terms of these mean 
fields; the inverse of the propagator S{k) can be read from 
the mean field Hamiltonian in the Nambu-Gor’kov bases 

IMllIS], 

i - M + /i 7 ° i 75 Z\ji?i 
- i75k \*Ri ^ - M - 
where the effective mass matrix has diagonal elements 

Mi = - 2Gs(7i + K\eijk\(7jak , (17) 

while the three diquark pairing amplitudes, 

A = -2d,7d, (18) 

and the effective chemical potential matrix, 

/i = /i — Gyn + fJ-sXs + /isAg -l- (J-qQ ■ 


(19) 

are color- and flavor-dependent. 

The thermodynamic potential for the quark part now 
reads 


[ 2 . 


ibare 


■Hbb = ^(qi'"qf ■ 


( 11 ) 


18 

= -2E 

i=i • 


d^k 

El +Tlnfl-be-l'^^l/^) 

J (27r)3 

[2 V J\ 


In principle we should include more general set of effec¬ 
tive interactions for other channels. We will keep only the 
vector part for simplicity. 

With this schematic Hamiltonian at hand, we will im¬ 
pose the constraints from charge and color neutrality, and 
/3-equilibrium conditions nil. For these conditions, we in¬ 
troduce three chemical potentials (Lagrange multipliers) 
and add the following terms to the Hamiltonian: 

^"Hconst = -f^giq^'Qq - llh) - Mal^Agg - /rsqUsq ,(12) 

where (/ig, ^ 3 , /xg) are chemical potentials for electric charge, 
the third and eighth components of the color matrices, re¬ 
spectively. (The other color charge densities automatically 
vanish for the condensates to be considered [551155] . 1 Q is 
the electric charge matrix Q = diag.(2/3, —1/3, —1/3), 
and = e,p) are electron and muon fields. 

The Hamiltonian will be treated within the mean field 
approximation. The mean fields for the chiral condensates 
and quark number density are 


3 

+ E] — 4iLtTi(T20’3 — 

i=l 


Gy 


nG(20) 


where the single particle contributions contain 18 inde¬ 
pendent eigenvalues, and A is the UV cutoff for the Dirac 
sea contribution. 

The ’’zero” of the renormalized thermodynamic poten¬ 
tial is set by 

T) = T) - = r = 0). (21) 

This choice of zero is demanded by the fact that the cos¬ 
mological constant be extremely small. 

The lepton contributions are standard. We add the 
renormalized (after subtraction of the vacuum energy) 
thermodynamic potential {I = e, /i) 

Qi = E / (1 + , (22) 


O'* = {q^qi), 




(13) 


Below we write (cti, 0 - 2 , 173 ) = {au,cFd,<ys) for later conve¬ 
nience. For the diquark mean fields, we write 


di = {q^Gi'y^Riq) , 


where 


(14) 


(15) 


with El = -f m/, and where we recall that the elec¬ 
tron chemical potential is = —fJ-Q. 

Writing the total thermodynamic potential as 17 = 
f2q + + 17^, the thermodynamic state of the system 

is determined minimizing the free energy with respect to 
the seven condensates {ai,di,n} under the neutrality con¬ 
ditions 

dG 

«Q,3,8 = - 

o/iQ,3,8 

which yields the “gap equations”, 
df2 dQ 
ddi 


= 0 , 


0 = - 




dai 


n = — 


dn 


(23) 


(24) 


Below we solve these self-consistent equations using the 
method outlined in [5911751. 


{Ri, R 2 , R 3 ) = (T 7 A 7 , tsAs, T 2 A 2 ) . 
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5.2 The impact of each interaction 

Our schematic model contains a number of parameters. 
For the NJL part, we fix the parameters {A,m, Gs, K) by 
the vacuum phenomenologj|2 To be specific, we use the 
set provided by Hatsuda and Kunihiro m- 

yl = 631.4 MeV, GsA^ = 3.67, KA^= 9.29, 
mu,d = 5.5 MeV, m, = 135.7 MeV. (25) 

This set gives Mu^d = 336 MeV and Ms = 528 MeV. For 
other choices of parameters, see 

The other parameters, H and Gy, are not constrained 
by the vacuum phenomenology. But from the viewpoint of 
the renormalization group, {H,Gv) are naturally the or¬ 
der of Gg. The standard choice is {H, Gy) = (0.75, 0.50)G's 
which is based on the Fierz transformation of the colored 
current-current interaction, ~ But this is just 

a tree level relation. The factor can easily be changed by 
the quantum corrections. Therefore we will consider the 
wider range for these parameters. 


5.2.1 The color magnetic interaction: H 


The standard choice of H is 0.75 Gs, while the choice H = 
Gg is called strong coupling. On the other hand, the choice 
H > Gg has been avoided, presumably because such large 
H yields the diquark condensates a t d ensity even lower 
than nuclear saturation density (Fig 13). In this case the 
2SC phase (Z\„d ^ 0, Ads = Agu = ^ appears first and 
then CFL phase {Aud, Ads, Agu ^ 0), but no normal quark 
matter phase. 

But the appearance of diquark condensates in very di¬ 
lute regime might be simply an artifact of non-confining 
models, for which the overlooked corrections are signif¬ 
icant at low density. Therefore, one should not take the 
problem literally and cannot use such problem to constrain 
the value of H which will be also used at high density. If 
we could add the confining effects to the NJL model, the 
colored diquarks would be combined with colored quarks 
to form baryons. In other words, in dilute regime diquarks 
would exist only inside of baryons, without showing them¬ 
selves in nuclear phenomenology. 

Moreover, there is a physics motivation to take H > 
Gg. Let us first note that in the NJL model the constituent 
quark mass for light quarks is ~ 336 MeV, larger than 
Mn/3 ~ 313 MeV. Thus the NJL pressure tends to appear 
at higher chemical potential than the hadronic pressure, 
and as the consequence these pressures look very differ¬ 
ent. But microscopically we should be able to describe 
hadronic matter in terms of quark degrees of freedom, 
so there must be some effects overlooked. One of the ef¬ 
fects to be implemented is attractive forces which reduce 
the average of quark energy to the smaller energy than 


® Strictly speaking, these parameters in medium can be 
smaller than the vacuum values. But as we will see, the neu¬ 
tron star constraints favor the values as large as the vacuum 


one. We will come back to this point in Sec 7.2 



[GeV] 



Fig. 13: Normalized baryon density ns/ris (top panel) and 
pressure P (bottom panel) as functions of quark chemical po¬ 
tential for several NJL parameter sets {Gy — 0 and sev¬ 
eral values for H) and APR. The only trustable region is 
shown with the bold lines {ns < 2ns for APR and ns > 5ns 
for NJL). At H > Gs, the baryon density and pressure in 
the NJL model begin to appear at lower chemical potential 
than Mjv/3 ~ 313 MeV where non-confining models are not 
trustable. 



0 0.5 1 1.5 2 

£ [GeV/fm^] 

Fig. 14: The pressure as a function of energy density. 
While H > Gs introduces the first order phase transition 
and the associated softening, at larger density the pairing 
effects eventually make the equations of state stiffer than 
those without the pairing. (See also discussions around 

Fig|) 
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Mat /3 as seen in the N — A splitting where the nucleon 
earns ~ 150 MeV energy reduction. We will use the color- 
magnetic interaction to describe this sort of effect, and 
our purpose is achieved by taking a large value, H > Gg- 

For a choice H > Gg, the quark matter pressure of 
non-confining models start to appear at /i < Mn/3 and 
can be larger than the hadronic pressure. The artificial 
excess of pressure should be suppressed by inclusion of 
the missing confining effects; packing quarks in the finite 
domain and the color flux tubes among quarks both cost 
energy. In this way of thinking, it is rather natural to have 
excess of low density pressure in the non-confining model. 
This picture motivates us to use H > Gg. The value of 
H is correlated with other parameters. If, for instance, 
the constituent quark mass appeared closer to (or smaller 
than) Mat/S, we could take smaller H\ if the quark mass 
appeared larger, we had to take larger H to achieve larger 
energy reduction. 

Finally we discuss the softening/stiffening induced by 
the color magnetic interactions. We have been accustomed 
to the examples in which the condensation phenomena 
soften the equation of state. This idea is not always true, 
however, as seen in Figj^ Certainly it is always true that 
the first order phase transitions soften equations of state, 
and most of model studies have emphasized only this as¬ 
pect. However, it is quite possible that the stiffness in the 
condensed phases grows fast and eventually exceeds that 
in the normal phase. This can be seen in FigfT^ where 
the diquark correlation induces overall shift of pressure 
curves toward lower chemical potential region, at least in 
the large ^ region. We have already discussed that such 
shift stiffens the equation of state (see discussions around 
Figj^. 

Therefore, the question relevant to neutron star phe¬ 
nomenology is at which density condensed phases become 
stiffer than the normal phase. For small values of H such 
reversal does not occur at density of interest, and this 
statement is consistent with the previous works. But for 
sufficiently large H, such reversal can occur at low den¬ 
sity as seen in Figj^ At large H{> 1.5), the system di¬ 
rectly goes to the condensed phase without experiencing 
the normal phase, and the condensed phase is stiffer than 
the normal phase from low to high densities. These dis¬ 
cussions show that whether condensed phases stiffen the 
equations of state or not is a quantitative issue. 


5.2.2 The repulsive density-density interaction: Gy 

The repulsive density-density interaction plays very im¬ 
portant roles to stiffen equations of state (Fig(^. It is 
intuitively clear that the repulsive force tries to prevent 
number density from increasing (Fig |16[ ); the repulsion in¬ 
creases the chemical potential whicn is necessary to in¬ 
crease the density and also pressure. The resulting P(/r) 
curves are slowly growing functions that are characteristic 
of stiff equations of state. 

Another important consequence of the vector interac¬ 
tion is that it smoothes out phase transitions [42] . Because 



0 0.2 0.4 0.6 0.8 1 1.2 1.4 

£ [GeV/fm^] 

Fig. 15: The pressure as a function of energy density. Large 
Gv leads to stiffer equations of state. 



[GeV] 



n [GeV] 


Fig. 16: Normalized baryon density nslus (top panel) and 
pressure P (bottom panel) as functions of quark chemical po¬ 
tential fj, for several NJL parameter sets {H = 0 and several 
values for Gy) and APR. 

baryon density changes slowly as chemical potential in¬ 
creases, the transition driven by baryon density also occurs 
slowly. As for the chiral transition at Gy = 0, the chiral 
restoration is triggered by the increase of baryon density, 
while the resulting reduction of the effective mass further 
accelerates the increase of baryon density. These two ef¬ 
fects drive each other and the changes occur very rapidly, 
resulting in the first order phase transition. On the other 
hand, increasing Gy brakes the increase of baryon den- 
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sity, so that the chiral restoration occurs only slowly. At 
large enough Gy, the vector interaction turns the first or¬ 
der phase transition into crossover; in general, introducing 
the repulsive density interactions shifts the chiral restora¬ 
tion line toward higher chemical potential region, and, as 
a result of the weaker strength of the phase transition, the 
location of the critical end point moves toward low tem¬ 
perature and eventually disappears at large Gy- In the 
standard NJL model, this occurs at Gy ~ OAGs- 

This feature suggests us an interesting interpretation 
on the lattice data for chiral restoration. The lattice Monte- 
Carlo calculations indicate that the curvature of the chiral 
restoration line at /r = 0 is unexpectedly flat m, and the 
chiral restoration line seems to be separated from the line 
characterizing the onset of baryon density. One of the sim¬ 
plest explanations on this behavior is to assume the exis¬ 
tence of repulsive density-density interactions. In NJL-like 
models, this effect can be expressed by the vector interac¬ 
tion. Seminal works [ZilTT] used this feature to constrain 
the value of Gy in the Polyakov-NJL (PNJL) model [58] . 
and they found that Gy jGg = 0.5—1.5 is reasonable range 
not to deviate from the curvature in the lattice data by 
^ 50%. (On the other hand, the authors in Pefs. [78117^ 
compared the quark number susceptibility in the effective 
models with the lattice result at T > Pc, and argued that 
Gy should be small enough not to suppress the suscepti¬ 
bility too much. In Sec |7.2[ we will discuss this case further 
in the context of gluon dynamics.) We will also see that 
a large Gy is necessary to get enough stiffness to sustain 
the 2 M 0 neutron stars. 

6 The 3-window equations of state 



Fig. 17: The 3-window pressures as functions of fj, for 
the parameter set {Gy,H) = (1.0,0.0)Gs (set 1) and 
(1.0,1.6)Gs (set 11). The set 1 gives unphysical interpola¬ 
tion: the curve contains the inflection point. At ub > 
the pressure with the set 11 almost coincides with the APR 
curve. The bold lines are used for APR at ns < 2ns and 
NJL at ub > 5ns, and their extrapolations are shown in 
thiner lines but with the same colors. 


6.1 The interpolated equations of state 

In this work we use a polynomial function to interpolate 
the APR and NJL equations of state. For the matching 
conditions at low and high densities, we demand P(^) to 
smoothly join the APR and NJL pressures up to the sec¬ 
ond derivatives; P, n, and the quark number susceptibil¬ 
ities are matched. Explicitly, first we define the chemical 
potentials at the boundaries as (ns = 3n)|^ 


In this section we construct 3-window equations of state 
by interpolating the APR and NJL equations of state. 
Then we calculate a neutron star mass to select out the 
equations of state which pass the 2Mq constraint. The M- 
R relation for non-rotating, spherically symmetric stars 
can be calculated through TOV equation, 


nB(fJ.<) = 2ns, nB{^J.>) = 5ns ■ (30) 

The polynomial function is 

5 

7^(M) = ^CfeM\ (31) 


dP _ GNis + P){M+4Trr^P) 
dr r(r — 2GnM.) 




(26) 


where we have six coefficients to satisfy three matching 
conditions at each boundary. 


, (32) 




d^P 


Q^Papr 


dip 




/ dr'A7:r'^e{r') . 

Jo 

(27) 




? 

M< 



dfjLi 



With the central baryon density n% (or central energy 
density Ec) as the boundary condition at r = 0 , we use 
the QCD equations of state and the TOV equations. The 


star radius R is determined as 


P(r = R-,n<B) = 0, 

(28) 

and the star mass is given by 


M{n‘B)=MiR-,n%). 

(29) 


For each equation of state we have a M-R curve as a 
function of n^. 


where j = 0 , 1 , 2 . 

As the form of the interpolating function suggests, here 
we are assuming the crossover behavior between hadronic 
and quark matter. But the 3-window construction itself 
does not necessarily demand such crossover; one can use 
the interpolating function which contains a kink which is 

® To check the sensitivity of our results to the boundary den¬ 
sities, we examined results by varying n< as 1.5 — and n> 
as 4 — 7no. These changes do not introduce any substantial ef¬ 
fects as far as the interpolation window is not taken to be too 
narrow. Also the effects of such changes can be easily absorbed 
by slightly changing the model parameters {Gy, H). 
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Q,/Qs = 0.5 QyQs = 0.8 G^Qs = 1.0 



^ [QeV] [QeV] n [QeV] 


Fig. 18: P v.s. for several vector couplings, GyjGs = 0.5,0.8,1.0. For each vector coupling, we take the value of H 
around which we can find physical interpolated equations of state. The curves with inflection points must be excluded. 
The bold lines are used for APR at ns < 2ns and for NJL at ub > 5ns. 


QyGs = 0.5 QyGs = 0.8 G,;Gs=1.0 



Fig. 19: ub/us v.s. fi. The curves with decreasing behaviors must be rejected. At high density, the curves with different H 
almost coincide, indicating that increasing H induces overall shift of pressure curves without changing the slopes much. 


characteristics of the first order phase transition, although 
its strength must be small enough to avoid significant soft¬ 
ening. In such case the final form of the 3-window equa¬ 
tions of state can be similar to those found in the conven¬ 
tional hybrid construction; the only difference in these two 
constructions are whether or not the extrapolated pres¬ 
sures are utilized explicitly. 

The polynomial function apparently allows many dif¬ 
ferent kinds of interpolated equations of state. Actually 
it is not so. Once we require the pressure curve to yield 
the physical speed of sound, the number of the candidates 
drastically decreases. 

Let us first consider the 3-window equation of state 
with large Gy, in order to pass the 2Mq constraint. While 
larger Gy certainly stiffens the equations of state, it also 
brings a trouble in sensible interpolations. As seen in Fig |17[ 
the NJL pressure with large Gy tends to appear at higher 
chemical potential region than the APR curve. As a con¬ 
sequence the NJL and APR curves cannot be connected 
without introducing the inflection points with which 
becomes negative. Thus we wish to increase Gy to pass 
the 2Mq constraint but the aforementioned unphysical 
behavior prevents us from doing so. 

The NJL pressure tends to appear at higher density be¬ 
cause the NJL constituent quark mass is Mg ~ 336 MeV, 
larger than the one-third of nucleon mass, ~ 313 MeV, 
that gives the threshold chemical potential of APR. As we 
saw in the last section, as we turn on the attractive pair¬ 
ing forces among quarks, the NJL pressure shows over¬ 
all shift toward the lower chemical potential region, ap¬ 


proaching the APR pressure. For the Gy = Gg case shown 
in Fig |17[ we increased the value of diquark coupling until 
we found the sensible interpolated pressure. We found that 
for Gy = Gg, the strong diquark coupling of ~ 1.6 G,, 
is necessary to avoid the inflection point. The resulting 3- 
window pressure almost coincides with the APR pressure, 
although they gradually depart each other as density goes 
beyond ^ 5ns. 

In Figj^we can see the set {Gv,H) = (1.0,1.6)Gs 
yields nonzero pressure at fj. < Mjq/Z. As we have empha¬ 
sized in several places, at low density, the extrapolated 
pressure of non-confining models may be (or should be) 
larger than the APR, since the confining effects are not in¬ 
cluded yet. As we will see, the equations of state passing 
the 2Mq constraint and speed of sound commonly show 
this excess pressure at low density. 

Now we move to more systematic studies in the param¬ 
eter space of (Gy, H). In Figs 18|19|20 we plot P, nsjns 
and Cg, respectively, as functions of /r. In Fig|^ we plot 
the neutron star mass as a function of the central number 
density, These results are given for Gy/Gg = 0.5,0.8 
and 1.0. As for the values of iJ, we searched the central 
value of H so as to give sensible interpolation at given 
Gy fGg, and then study the slight variations from it. Sev¬ 
eral remarks are in order: 

(i) For sensible interpolations, the appropriate values 
of H becomes larger as we increase Gy. For GyjGg = 
0.5,0.8, and 1.0, the appropriate choices of H/Gg are 1.4, 
1.5, and 1.6, respectively (we have not tried to fine-tune 
the values of H more than the accuracy of 0.1). To check 
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Q^Gs = 0.5 QyQs = 0.8 QyQs = 1.0 



Fig. 20: The (square of) speed of sound Cg as a function of fi. The negative exist but are not plotted. For each Gv, only one 
of curves (bold, red) has the physical speed of sound. (Strictly speaking, at Gv = Gs the value of Cg is slightly larger than 1 
and we need to take H slightly smaller than 1.6 to pass the causality constraint. This change does not affect our arguments.) 


G,/Qs = 0.5 = G,/Gs=1.0 



Fig. 21: The neutron star mass M/Mq v.s. ns/ris for the parameter set discussed in Fig 18 Only bold (red) lines are physically 
acceptable. To pass the 2Mq constraint, the vector coupling needs to be sufficiently large, Gv jGg > 0.5. 


whether the curves are physical or not, we examine the 
behaviors of number density and speed of sound in FigsfT^ 
and[20l 

(ii) The resulting pressure curves are similar to the 
APR pressure, indicating that the boundary conditions 
up to the second derivatives are considerably tight. The 
matching domain becomes wider for larger Gy, beyond 
such domain the NJL pressure gets softer than APR. It 
seems that Gy/G^ = 0.5 is a bit too small to pass the 
2 M 0 constraint, while Gy/Gs = 0.8 and 1.0 allows us to 
pass the constraint. 

(iii) The conditions on the speed of sound appear to 
be quite tight. Because the APR is relatively soft at = 
1 — 2 ns, the stiff equations of state at high density side nec¬ 
essarily requires the rapid growth in P{e). With careless 
interpolation one easily hnds the violation of the causality. 

(iv) As seen from Fig|M] a neutron star with the mass 
M/Mq > 0.7 has the central density of ub > 2ns, and 
requires information beyond purely nuclear matter. The 
equations of state at ub > 4ns carry vital information 
about whether they can pass the 2Mq constraint. 


6.2 The mass-radius relations 

For the M-R relations, we pick up only physical 3-window 
equations of state from Figj^ the parameter sets are 
(Gy/Gs,iT/Gs) = (0.5,1.4), (0.8,1.5), and (1.0,1.6). The 
results are shown in Fig {^together with the lines of 1.97± 



R [Km] 

Fig. 22: The M-R relation. From the equations of state 
displayed in Fig{^ we picked out only physically accept¬ 
able ones. For canonical neutron stars with M ~ IAMq, 
the radii are about R ~ 11.2 —11.5 km, mainly determined 
by APR for ns ^ 2ns. 


O.O 4 M 0 for J1614-2230 and of the canonical mass, IAMq. 
As we already discussed, larger Gy leads to the larger 
neutron star mass. As for the neutron star radii, R ~ 
11.2 — 11.5 km at the canonical mass. 

The shape of a M-R relation deduced from our 3- 
window result is similar to that from the APR, as ex¬ 
pected from the behaviors of pressure. But we should no- 
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tice that the degrees of freedom in the equations of state 
are quite different; while the APR contains only nucle¬ 
onic degrees of freedom to very high density, the 3-window 
equations of state manifestly take into account the trou¬ 
blesome strangeness degrees of freedom. This difference 
should be crucial since the 2Mq neutron stars need infor¬ 
mation around ^ Srig where strangeness is supposed to 
exist. 

The consequences of our 3-window and conventional 
hybrid equations of state can be distinguished by the neu¬ 
tron star radii. In conventional hybrid construction, we 
tend to choose stiff hadronic equations of state at low den¬ 
sity that enables stiff quark matter pressure to intersect 
with the hadronic one from below. With the stiff hadronic 
equations of state, the neutron star radii tends to be large, 
i? = 13 — 16 km. In contrast, in our 3-window construc¬ 
tion we use soft hadronic equations of state at low density, 
and, at the same time, stiff quark matter equations of state 
at high density; the former is responsible for small radii 
of ~ 11 km while the latter is used to yield the mass of 
~ 2Mq. 

7 Discussions 

So far we have examined the interpolated pressure but 
have not discussed the physics in it. In this section we will 
infer the physics in the interpolated region. In the follow¬ 
ing we will discuss the hyperon puzzle and the nature of 
gluon dynamics at density relevant to neutron stars. 


7.1 On the hyperon problem 

In purely hadronic descriptions with hyperons, the many- 
body approaches for potential models such as Hartree- 
Fock [T7], Brueckner-Hartree-Fock [TS] give equations of 
state which encounter the hyperon softening around ns = 
2 — 3ns. The resulting maximal star mass is about ~ 
I.6M0, smaller than the 2Mq constraint. The results in¬ 
dicate that for the NN and YN potentials {N and Y are 
nucleons and hyperons, respectively), not only repulsive 
two-body forces but also more-body repulsion are neces¬ 
sary to prevent matter from softening. Recent Quantum 
Monte Carlo calculations showed the big impact of the 3- 
body forces which have yet been determined empirically 

HO]. 

The mechanisms proposed in the context of hadronic 
descriptions typically utilize the higher order effects such 
as many-body forces to delay the appearance of hyper¬ 
ons. But whenever the higher order effects have too much 
importance, there should arise the questions about the 
convergence. Then the relevant question to be answered 
is how they behave after summing up the higher order 
corrections. 

To guess how the results change toward high density, 
let us consider the hyperon problem from the microscopic 
viewpoint. In quark model descriptions, the strangeness 
begins to appear when quark chemical potential reaches 
Ms ~ 0.5 GeV if we ignore the chiral restoration. In terms 


of baryon chemical potential, it amounts to the order of 

~ 3Ms ~ 1.5 GeV, considerably higher than the A 
or A-masses, ^ + Ms ~ 1.1 GeV. Thus the 

strangeness appears at relatively higher chemical potential 
in the quark description than in the hadronic one. 

Let us consider why this is the case. Let us first use 
purely hadronic descriptions. For simplicity, we consider 
pure neutron matter with A hyperons and ignore the in¬ 
teractions for the moment (therefore the numerics below 
should not be considered too seriously). In purely hadronic 
descriptions, A at rest appears through a weak decay, 
n{udd) —^ A(uds), when the neutron Fermi energy EnP = 
m'^ + reaches the A-mass. The corresponding neu¬ 
tron Fermi momentum is ~ 0.6 GeV ^ Aqcd- Treating 
baryons as if elementary particles, the onset baryon chem¬ 
ical potential for A is /is = vriA- 

The above description, however, ignores important as¬ 
pects of the matter made of composite particles. At baryon 
chemical potential where the neutron Fermi momentum 
becomes the order of Aqcd, the phase space for u and 
d-quarks is occupied to quark momenta of ^ Aqcd- In 
such situation, A baryons, composite particles made of 
u, d and s-quarks, cannot be free from the Pauli blocking 
effects acting on u, d quarks. For instance, A baryons at 
rest are dominantly described by quarks at low momenta, 
and their emergence is obstructed by already existing u, d 
quarks in the neutron Fermi sea. While the A-baryons at 
sufficiently large momenta can escape from these Pauli 
blocking effects, such A is too energetic and cannot be 
generated through the neutron weak decays. One can cre¬ 
ate A at rest by replacing d-quark at low momenta with 
s-quark, but it increases the energy of the system by ~ 
Ms — Md, thus cannot occur in the /3-equilibrium. The 
only strange baryon that is perfectly free from the quark 
Pauli blocking is f2~{sss), but it appears only when fip 
reaches ma 1.6 GeV ^ 3Ms. The last case roughly cor¬ 
responds to the onset baryon chemical potential in quark 
matter descriptions. 

To be more concrete, consider a baryon in a constituent 
quark model, for example, of Isgur and Karl [5^. The 
spatial part of the ground state baryon wavefunction at 
momentum K can be expressed by quark wavefunctions 
as (P = pi -f P2 +P3) 


'I'T'^\puP2,P^) = 


(27r)3/2 


5{K-P) 



3 9x2 \t'p ' t'K.f 


,(33) 


where V is the volume of the system, A ^ characterizes the 
size of wavepacket, and Pp and are Jacobi’s momenta. 


Pp = 


Pi -P2 

V2 


Pk 


P 1 +P 2 - 2p3 

V6 


(34) 


The amplitude is large when pi ~ P2 P3 ~ P/3. Thus if 
a baryon has the momentum IT ^ 0, the quarks inside of 
it occupy the phase space at low momenta. At low baryon 
density, low momentum states of quarks are only partially 
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filled, so baryons do not feel the Pauli blocking effect at 
quark level. As the baryon Fermi sea develops, however, 
the quark states at low momenta are largely occupied, and 
a baryon having those states as a part is obstructed; to 
introduce another baryon into the system, we have to keep 
the baryon’s momentum large enough to avoid conflicts 
with the quark Pauli blocking. 

To quantitatively understand whether the quark Pauli 
blocking becomes significant at density relevant for neu¬ 
tron stars, more details must be worked out by includ¬ 
ing all quantum numbers (color, spin, flavor) and interac¬ 
tions for baryons. But the above arguments at least illus¬ 
trate that the hyperons in microscopic descriptions do not 
emerge as easily as in purely hadronic models. 

Within purely hadronic descriptions, we expect that 
the aforementioned quark Pauli blocking effects can be 
mimicked by the effective repulsive forces among baryons. 
This might justify phenomenological introduction of re¬ 
pulsive YN forces. But even so, such treatments work only 
partially because our ability to handle many-body forces 
is limited. In contrast, the arguments at quark level have 
advantages at large density where a hadronic description 
begins to reveal its weakness. 


7.2 On the nature of gluon dynamics 

So far our model descriptions have not manifestly taken 
into account the gluonic degrees of freedom. Some of the 
gluons may be integrated out into the parameters such as 
Gs, Gy, and H in the NJL model. Those parameters are in 
principle density dependent and should become smaller at 
larger density as the gluons become weakly coupled due 
to medium screening effects. In this respect, the impact 
of the 2Mq constraint is substantial, since it requires the 
NJL parameters to be as large as Gs in vacuum. As we 
saw in the last section, GyjGs > 0.5 and H/Gs > I.O 
seems necessary for the region of ^ 5 — lOn,, that 
corresponds to ^ ^ 0.5 — 0.7 GeV. This suggests that at 
density relevant to neutron stars the gluon sector remains 
strongly non-perturbative as in the QCD vacuuirp^ 

Actually this picture is consistent with another im¬ 
plicit assumption in our model treatments on the QCD 
vacuum energy. 

It has been argued that the full QCD vacuum has 
smaller energy than the perturbative vacuunj^ In the 
NJL model, a part of vacuum energy has been automati¬ 
cally taken into account through the Dirac sea contribu¬ 
tions associated with the chiral symmetry breaking, (to^: 
current quark mass, constituent quark mass) 

B, = Gginig) - n.iMg) ^ > 0, (35) 

This point of view was more manifestly examined in 
Ref.|81|. where the authors discussed how the medium depen¬ 
dent Gv should behave as a function of /r. 

The terminology such as perturbative and non- 
perturbative vacua might sound awkward as mentioned 
in Ref. [S2], but we will keep using it as long as it does not 
affect our arguments substantialy. 


which we call the quark bag constant. On the other hand, 
this is not the full contribution to the vacuum energy, since 
non-confining models like NJL do not take into account 
the long range components of gluons; those gluons, upon 
integration, would leave non-local interactions in the ef¬ 
fective models for quarks, but such interactions are not 
present in the NJL modef^ Therefore besides the Dirac 
sea contributions, there should be additional energy re¬ 
duction in the QCD vacuum that is purely caused by the 
gluons in the infrared. We call it the gluonic bag constant, 
Bg > 0, which are related to the gluon condensate. These 
positive energy contributions should appear whenever the 
non-perturbative matter changes into perturbative mat¬ 
ter; when applying perturbative calculations we must add 
the positive constant term (bag constant) to the pertur¬ 
bative energy density because such calculations perform 
an expansion around the trivial vacuum, not around the 
true QCD vacuum (whose energy density is renormalized 
to zero at = T = 0). 

Now imagine that the gluons get screened at ns ~ 
2 — lOus, and gluons become weakly coupled. Then the 
size of the gluon condensate is reduced and the gluonic 
bag constant must be added to (subtracted from) our NJL 
energy (pressure). But the addition of such constant can¬ 
not be as large as ~ , because it would substantially 

soften the equations of state in which the NJL pressure is 
only about 0.2 — 0.5GeV/fm^ ~ 2/1^^^ at ub ~ 5ns. 
Clearly it is no longer possible to pass the 2Mq con¬ 
straint. Strictly speaking, even in such situation one can 
still achieve the necessary stiffness by taking Gy and H 
much larger than the values we took, but the appearance 
of the large bag constant and substantially large NJL pa¬ 
rameters at finite density seem to be incompatible assump¬ 
tions. Rather, it is much more consistent to take Gy and 
H as the values we took, and then also keep the gluonic 
bag constant to the vacuum value which is much smaller 
than ~ Aqqq. 

In this view on the medium NJL couplings, our NJL 
parameters (Gg, AT, Gy, H) should decrease when the pres¬ 
sure curves approach the pQCD domain. Otherwise the 
NJL pressure undershoots the pQCD results and the speed 
of sound is larger. This requirement is again consistent 
with the idea that the NJL couplings originate from non- 
perturbative gluons and at very large density should dis¬ 
appear together with such gluons. 


Finally we comment on the estimate on Gy which 
was briefly mentioned in Sec 5.2.2 In Ref. the au¬ 
thors performed the PNJL studies for quark number sus¬ 
ceptibilities at T > Tc, and showed that those quanti¬ 
ties are suppressed too much in the presence of Gy. On 
the other hand, they explained the flatness of the chiral 
restoration curvature using hadronic models with repul¬ 
sive density-density interactions for T < T^. But it should 
be able to express the hadron-hadron interactions in terms 
of quark degrees of freedom, so we will not distinguish 


This two scale argument is analogous to those of Manohar 
and Georgi [83] , and Weinberg [84| . 
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hadron-hadron and quark-quark interactions rigidlj[^ in 
the latter description the repulsive density-density inter¬ 
action at T < Tc are effectively mimicked by large Gy. 
Thus as explained in Ref. m, our interpretation is that 
Gy is large aX T < and small at T > Tc. Then we next 
notice that one of the major characteristic difference be¬ 
tween T <Tc and T > T^is the gluon dynamics. Since we 
get the NJL parameters after integrating out gluons, the 
changes of the NJL parameters such as Gy around T ^ Tc 
are rather natural. We also point out that the changes of 
the NJL parameters driven by gluons should be universal 
for all different flavors, and such picture seems to be con¬ 
sistent with the lattice studies on the off-diagonal flavor 
correlators m- 


8 Summary 

In this paper we discuss the 3-window construction of 
QCD equations of state, and explain why this construc¬ 
tion can express the physics that are largely unexplored in 
the conventional hybrid or self-bound quark matter equa¬ 
tions of state. In the 3-window picture, we have purely 
hadronic matter, percolated quark matter, and matter in¬ 
termediate between these two. Such matter is presumably 
expressed as either hadronic matter about to percolate 
or quark matter about to be confined. The range of such 
matter is ~ 2 — bUg or ^ ~ 0.35 — 0.6 GeV. 

We examine the significance of the 2Mq constraint 
which rules out soft high density equations of state. The 
constraint is particularly severe when the purely hadronic 
equations of state are soft: in general it is more difficult 
to construct equations of state that are soft at low density 
and stiff at high density, because the combination of them 
tends to violate the causality constraint. But according to 
recent results from many-body calculations |29j and im¬ 
plications of rather small neutron star radii [Toiiiiiiiaiia 
[Hiin], it seems necessary to construct such soft-stiff equa¬ 
tions of state. We have explained that the 3-window con¬ 
struction has more chances to provide the soft-stiff equa¬ 
tions of state than the conventional hybrid one with the 
first order phase transition, although the final form of the 
equations of state can be similar if the strength of the 
phase transition is very weak. 

The 3-window picture also changes the guidelines for 
practical modeling. Within the NJL model, we found that 
Gy > 0.5Gs and H > Gg are necessary to pass the 2Mq 
constraint with the physical speed of sound, while quark 
model pressure apparently exceeds the hadronic pressure 
at low /r. In contrast to the conventional hybrid construc¬ 
tions, we accept such behaviors by expecting that at low 
density the confining effects largely suppress the artifi¬ 
cial excess of pressure in the non-confining models. With 
such expectation, we even think such excess pressure to be 

For instance, the channel dependence of the nucleon- 
nncleon interactions can be understood by just summing 
up quark exchange contributions. They can be organized as 
1/Nc ~ 1/10 expansion. See Refs. |85ll86 | and the references 
therein. 


quite natural; if the pressure of non-confining models has 
less pressure than the hadronic one, we must imagine that 
the confining effects enhance the pressure to approach the 
full result — but this is unlikely. 

The values of Gy and H used for the percolated quark 
matter at > Srig are as large as Gg in the vacuum. Since 
these NJL parameters should be obtained after integrat¬ 
ing out gluons, the large NJL couplings in the percolated 
domain indicate that the nature of gluons are not changed 
appreciably by the medium effects to ~ lOn^. This im¬ 
plies that at ns ~ 5 — lOris we have the percolated quark 
matter with non-perturbative gluons which would try to 
keep the system locally color singlet. The useful regime to 
study such matter is large Nc^ as discussed by McLerran 
and Pisarski [45] . 
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